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ABSTRACT 

We study, by means of dedicated simulations of massive black hole build-up, the 
possibility to constraint the existence and nature of the AGN population at z > 6 
with available and planned X-ray and near infrared space telescopes. We find that 
X-ray deep-field observations can set important constraints to the faint-end of the 
AGN luminosity function at very high redshift. Planned X-ray telescopes should be 
able to detect AGN hosting black holes with masses down to > 10^ Mq (i.e., X- 
ray luminosities in excess of 10^^ erg s~^), and can constrain the evolution of the 
population of massive black hole at early times (6 < z < 10). We find that this 
population of AGN should contribute substantially (~ 25%) to the unresolved fraction 
of the cosmic X-ray background in the 0.5-10 keV range, and that a significant fraction 
(~ 3 — 4%) of the total background intensity would remain unaccounted even after 
future X-ray observations. As byproduct, we compute the expected UV background 
from AGN at z > 6 and we discuss the possible role of AGN in the reionization of the 
Universe at these early epochs, showing that AGN alone can provide enough ionizing 
photons only in the (improbable) case of an almost completely homogeneous inter- 
galactic medium. Finally, we show that super-Eddington accretion, suggested by the 
observed QSOs at z ~ 6, must be a very rare event, confined to black holes living in 
the highest density peaks. 
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1 INTRODUCTION 

The formation of black hole (BH) seeds and their evolu- 
tion have been the subject of several theoretical investiga- 
tions. The "flow-chart" presented Rees (1978) still stands as 
a guideline for the possible paths leading to formation of BH 
seeds in the center of galactic structures. One possibility is 
the direct formation of a massive BH (MBH) from a collaps- 
ing gas cloud (Haehnelt & Rees 1993; Loeb & Rasio 1994; 
Eisenstein & Loeb 1995; Bromm & Loeb 2003; Koushiappas, 
Bullock & Dekel 2004; Begelman, Volonteri & Rees 2006). 
The gas can condenses to form a central massive object. The 
mass of the seeds predicted by different models vary, but typ- 
ically are in the range Mbh ~ iC — 10® Mq . Alternatively, 
the seeds of BHs can be associated with the remnants of the 
first generation of stars, formed out of zero metallicity gas. 
The first stars are believed to form at z ~ 20 in halos which 
represent high-a peaks of the primordial density fie Id. The 
absence of metals might lead to a very top-heavy initial 
stellar mass function, and in particular to the production 
of very massive stars with masses > lOOAf© (Carr, Bond, & 



Arnett 1984). If very massive stars form above 260 Mq, they 
would rapidly collapse to BHs with little mass loss (Fryer, 
Woosley, & Heger 2001), i.e., leaving behind seed BHs with 
masses Mbh ~ 10^ - 10^ Mq (Madau & Rees 2001). The 
subsequent growth of BHs from these initial seeds has been 
investigated in the hierarchical framework, typically associ- 
ating episodes of accretion to galaxy mergers (Haiman 2004, 
Yoo & Miralda-Escude 2004, Shapiro 2005; Volonteri & Rees 
2005, 2006; Lapi et al. 2006). 

Only a few observational constraints are currently avail- 
able. Many are limited to the brightest sources, probing only 
the very bright end of the luminosity function of quasars at 
z ~ 6 and therefore only the upper end of the massive black 
hole mass function (Fan et al. 2001). The observation of very 
luminous quasars, powered by billion solar masses BHs, at 
z « 6 in the SDSS survey (Fan et al. 2001), implies that a 
substantial population of smaller accreting BHs must exist 
at earlier times. Deeper but limited constraints have been 
provided for more typical AGN, with ~ 2 orders of mag- 
nitude smaller MBHs, using deep X-ray observations (e.g., 
Alexander et al. 2001; Barger et al. 2003b; Koekemoer et 
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al. 2004; Mainieri et al. 2005). However, these constraints 
lack the detail and depth required to understand and de- 
termine the global evolution of BH seeds and of the aver- 
age MBH population at very high redshifts. Accreting BHs 
arc observed as X-ray emitters up to 2 < 6 (Barger et al. 
2003, Steffen et al. 2006), and there are no reasons to believe 
that higher redshift AGN would be any different. The detec- 
tion of very high redshift AGN is one of the goals of future 
space missions. In particular. X-ray telescopes, such as the 
planned XEUS and Constellation-X, arc expected to detect 
these sources even at « > 6. Moreover, X-ray observatories 
can identify even heavily obscured AGN activity, due to the 
penetrating nature of hard X-rays. 

Additional, albeit less direct, limits on the early MBH 
population can already be placed by the requirement that 
the cumulative emission from the predicted high redshift 
sources does not saturate the observed unresolved X-ray 
background (Dijkstra, Haiman & Locb 2004; Salvaterra, 
Haardt, Ferrara 2005). We show that future space-borne 
missions will be able to constrain different proposed models 
for the accretion history of MBHs at early times. We focus 
here on a specific model for MBH formation, which traces 
MBH seeds to the first generation of metal-free stars. Wo 
then compare how different accretion histories reflect onto 
the detectability of X-ray sources. 

The mass growth of the most massive BHs at high red- 
shift must proceed very efficiently to explain the luminosity 
function of luminous quasars at z « 6 in the Sloan Dig- 
ital Sky Survey (SDSS, e.g. Fan et al. 2001). Volonteri & 
Rees (2005) explore the conditions which allow a sufficient 
growth of MBHs by 2 = 6 under the assumption that accre- 
tion is triggered by major mergers. At such high redshift an 
investigation of evolution of the MBH population must take 
into account the dynamical evolution of MBHs. Black hole 
mergers, in fact, can give a net negative contribution to the 
early black hole growth, due to gravitational effects which 
can kick BHs out of their host halos (gravitational recoil, 
due to to the non-zero net linear momentum carried away 
by gravitational waves). Adopting recent estimates for the 
recoil velocity (Baker et al. 2006), Volonteri & Rees (2005) 
find that if accretion is always limited by the Eddington rate 
via a thin disc, the mgiximum radiative efficiency allowed to 
reproduce the LF at 2 = 6 is €max = 0.12 (correspond- 
ing to an upper limit to the MBH spin parameter of 0.8). 
If, instead, high-redshift MBHs can accrete at super-critical 
rate during an early phase (Volonteri & Rees 2005, Begel- 
man, Volonteri & Rees 2006), then reproducing the observed 
MBH mass values, is not an issue. The constraints from the 
LF at 2 = 6 are still very weak, so either a model with a 
low emax or a model with supercritical accretion cannot be 
ruled out based on these results only. 

In a previous paper we have investigated how cur- 
rent observations can constrain the late accretion history 
of MBHs, at z < 4 (Volonteri, Salvaterra & Haardt 2006). 
In this paper we study by means of a detailed model of MBH 
assembly, the detectability of very high redshift AGN both 
at X-ray and near-IR wavelengths by future missions. We 
discuss the contribution of these objects to the unresolved 
X-ray background (in particular in the 0.5 2 keV band). 
Finally, we discuss the role of an early population of AGN 
in the reionization of the Universe and their contribution 
to the UV background at high redshift. The paper is orga^ 



nized as follows. In Sect. 2 we describe briefly the merger 
tree model for the formation and evolution of MBH in the 
early Universe. Sect. 3 presents the basic equations we need 
in the number counts and background calculation and, in 
Sect. 4, we discuss the adopted spectrum of AGN. In Sect. 5 
we present our results and discuss the role of future space 
mission in the observations of high redshift AGN. In Sect. 6 
we compare our results to those of a model in which the 
early evolution of MBH seeds is characterized by a phase of 
super-critical accretion. Finally, in Sect. 7, we briefly sum- 
marize our results. 

Throughout the paper we use the AB magnitude 
system* and the standard ACDM cosmology (Spergel et al. 
2003). 



2 HIGH REDSHIFT MASSIVE BLACK HOLE 
EVOLUTION 

The main features of a plausible scenario for the hierarchical 
assembly, growth, and dynamics of MBHs in a ACDM cos- 
mology have been discussed in Volonteri, Haardt & Madau 
(2003). Dark matter halos and their associated galaxies un- 
dergo many mergers as mass is assembled from high redshift 
to the present. The halo merger history is tracked back- 
wards in time with a Monte Carlo algorithm based on the 
extended Press-Schechter formalism. "Seed" holes are as- 
sumed to form with intermediate masses in the rare high 
V — a peaks collapsing at 2 = 20 — 25 (Madau & Rees 2001) 
as end-product of the very first generation of stars. As our 
fiducial model we take u = 4: &t z = 24, which ensures 
that galaxies today hosted in halos with mass larger than 
10^^ Mq are seeded with a MBH. The assumed threshold al- 
lows efficient formation of MBHs in the range of halo masses 
effectively probed by dynamical studies of MBH hosts in the 
local universe. 

As a reference, we adopt here a conservative model as- 
suming Eddington-limited accretion which is able to repro- 
duce the bright end of the optical LF, as traced by obser- 
vations in the SDSS (e.g. Fan et al 2001). We then dis- 
cuss in Section 6 a simple model which considers super- 
Eddington accretion rates for high-redshift MBHs and one 
in which more massive seeds form late, as in Koushiappas et 
al. (2004), and evolve through Eddington-limited accretion. 

The fraction of active BHs (i.e., AGN) at early times is 
large in our models, reaching almost unity for the most mas- 
sive black holes (Mbh > 5 x 10^ Mq) at 2 ~ 5. The fraction 
of AGN with respect to the whole BH population decreases 
by at least one order of mag nitudo at Mbh ~ 10^ - lO'' Mq 
and further at lower masses. Under our model assumptions, 
detection of active sources can therefore provide robust con- 
straints on the whole BH population. 



* AB magnitudes are defined as AB= — 2.5 log(F,^())— 48.6, where 
is tlie spectral energy density within a given passband in units 
of erg cm~^ Hz~^ 
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3 BASIC EQUATIONS 

The number of sources observed per unit solid angle at red- 
shift Zo in the flux range F^^ and F,^g + dF„g at frequency 

vo is 



dN 



dndF^ 



-{F^o,zo) 



(1) 



where dVc/dzdfl is the comoving volume element per unit 

redshift per unit solid angle, and nc{z, F^^) is the comoving 
density of sources at redshift z with observed flux in the 
range , F^^ +dF,^g]. The specific fiux of a source observed 
at Zo is given by 



L,.(MBH)e 



(2) 



where L,^{Mbh) is the specific luminosity of the source (in 
units of erg s~^ Hz~^) averaged over the source lifetime, 
which is assumed to be only a function of the mass of the 
central BH. In the above cq. (2), u = h'o{l + z)/{l + zq), 
dL{z) is the luminosity distance, and Teff is the effective op- 
tical depth of the intergalactic medium (IGM) at fo between 
redshifts zo and z. Shortward of the observed Lyman— a, the 
emergent spectrum is strongly modified by IGM absorption 
(see Sect. 2.2 of Salvatcrra & Ferrara 2003 for a full descrip- 
tion of the IGM modeling). For sources at a > 6, as those 
we are considering in this paper, IGM absorption in the ob- 
served 0.5-10 keV band (corresponding, for sources at 2: > 6 
as those considered here, to rest-frame energy > 3.5 keV) 
can be safely neglected. 

The radiation background Jvg{zo) observed at redshift 
Zo at frequency vo, is 



(1 + zo) 



77 



(3) 



where ev{z) is the comoving specific emissivity, and dl/dz is 
the proper line element. The source term tv is given by 



e.{t) 



= [ dMsH [ L4t-t',MBH) dt' 
J Jq dt'dMBH 



dMBH rL,.(MBH) 



dtdMBB ' 



(4) 



The second approximated equality holds once we consider 
the source light curve averaged over the typical source life- 
time T, and assuming the source formation rate per unit 
mass as constant over such timescale. 



4 AGN SPECTRUM 

The physical characterization of the source is encoded in 
its spectral energy distribution (SED), L,^. We model the 
UV part of the SED of unabsorbed AGN (log(iVH/cm-2) < 
22, referred as type I) as a multicolor disk black body 
(Shakura & Sunyacv 1973). Assuming Eddington limited 
accretion, the maximum disk temperature is fcTmax ~ 
1 keV{MBH /Mq)~^^'^. The characteristic multicolor disk 
spectrum is broadly peaked at Epea.k ^ SkTmax, and fol- 
lows a power law with Li, cx u^^^ for energies hu < Speak, 
and exponentially rolls off for E > E'peak- In the X-ray, 
the spectrum can be described by a power-law with pho- 
ton index F = 1.9 {f^ oc i'"'^), and an exponential cutoff at 



Ec = 500 keV (Marconi ct al. 2004). The averaged X-ray 
SED of absorbed AGN (type II) is described by the same 

type I spectrum for iJ > 30 keV, and, in the range 0.5 30 

keV, by a power-law (continuously matched) with photon 
index F — 0.2, obtained convolving the type I spectrum 
with a lognormal distribution of absorption column density 
centered at Iog(iV/f/cm"^) = 24 (Sazonov et al. 2004). UV 
emission from type II AGN is assumed to be negligible. We 
further assume a type I/type II ratio of 1/4, independently 
of redshift and luminosity, since detailed modelling of in- 
trinsic absorption at « > 6 are not currently available. Note 
that the observed X-rays correspond, for sources at « > 6, 
to rest frame energies between 3.5 and 70 kcV, where the 
emission properties of type I and type II AGN are thought 
to be similar (the rare Compton thick sources are not con- 
sidered here). 

The X-ray emission of type I AGN is normalized to 

the optical, adopting an optical-to-X rays energy index 
aox = 0.126 log(L25()o) + O-Olz - 2.311 (sec Eq. 5 of Stef- 
fen et al. 2006), where L2500 is the monochromatic lumi- 
nosity at A = 2500 A(rest frame). By definition, the X- 
ray luminosity at the rest frame energy of 2 keV, L2, is 
L2 = i25()()(i^2/i^2r)()()) The scaling of aox with redshift 

and luminosity has been obtained by Steffen et al. (2006) 
combining data from the SDSS, COMBO-17 and Chandra 
surveys. The final sample consists of 333 AGN extending out 
to a ~ 6 spanning five decades in UV luminosity and four 
decades in X-ray luminosity. Since only a mild dependence 
on redshift is found, we extrapolate this result also to 2; > 6. 
We further assume Eddington limited accretion (see Section 
2). 



5 RESULTS 

5.1 X-ray number counts 

We compute the soft (0.5-2 keV) and hard (2-10 keV) X- 
ray number counts from MBHs shining at z > 6 predicted 
by our model of MBH assembly and evolution. Results are 
shown in Fig. 1, and arc compared to available observa- 
tional data (Moretti et al. 2003). The bow-tie indicates re- 
sults from the fluctuation analysis of the Chandra deep field 
(Miyaji & Griffiths 2002). Different fines refer to different 
ranges in BH masses. At the fiux limit of current surveys, 
log(S') > -16.6(-15.8)crgs^^cm"^ in the soft (hard) band 
(Alexander et al. 2003), the contribution to the logN/logS 
from sources at 2 > 6 is ~ 8% in the soft X-rays, < 1% 
in the hard band. In the 0.5-8 kcV band, we expect ~ 2 
source at ^ > 6 in the Chandra Deep Field North with 
fluxes exceeding 3 x 10~^^ erg cm~^, to be compared 
to an upper limit of 7 sources with extreme X-ray /optical 
flux ratios (EXO, Koekemoer et al. 2004), that are candidate 
high redshift AGN. 

The optical identification of these objects is problem- 
atic, owing to absorption (both internal and by the IGM) 
and to their very low optical flux. Barger et al. (2003a, 
2003b) searched for optical counterparts (at 5a confidence) 
with z band magnitude zaso < 25.2 of AGN of the Chandra 
Deep Field North exposure. Apart from a source at z = 5.19 
with 2860 = 23.9, no other 2 > 5 candidate was identified. 
This is consistent with our results, where the majority of the 
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Figure 1. Predicted logN/logS in the soft (upper panel) and 
hard {lower panel) X-ray bands for sources at z > 6. Different 
lines refer to different BH mass ranges: all masses (solid line), 
Mbh = lO"* - 10^ Mq (dotted line), Mbh = 10^ - 10^ Mq 
(dashed line), Mbh = 10^ — lO'^ Mq (long dashed line), and 
Mbh > 10'' Mq (dot-dashed line). Data points are the measured 
source counts obtained by Chandra (Morctti ct al. 2003) and the 
bow-tie indicates results from the fluctuation analysis (Miyaji & 
Griffiths 2002). 

very high redshift sources detected in the deepest Chandra 
observations should have a z-band magnitude fainter than 
27. Barger et al. (2003b) used the lack of optical identifi- 
cations to derive limits on the number of objects at z > 5. 
They found that in a field of view (f.o.v.) corresponding to 6 
arcniin radius circle, only ~ 6 sources with fluxes exceeding 
2 X lO"^*^ erg s"^ cm"^ in the 0.4-6 keV band should lie at 
« > 5. Our model is consistent with this limit, predicting 
that only ~ 1 sources at 2: > 6 should be found in the same 
field of view. 

At fainter fluxes, a signiflcant fraction of the sources 
identified in the fluctuation analysis of the deepest Chan- 
dra data might be AGN at z ~ 6. Our model predicts that 
high redshift AGN contribute to the logN/logS at a level 
~ 13 — 60% in the soft X-ray band, < 6% in the hard X-ray 
band at the flux limit of the fluctuation analysis (~ 8, 000 
sources deg~^ at the flux limit logS = —17.2 for the soft 
band, and ~ 300 sources deg~^ at the flux limit logS = — 16 
for the hard band). Direct observations of such sources is 
among the main goals of the next generation of X-ray tele- 
scopes (e.g. XEUS^ and Constellation-X^). The XEUSraSs- 
sion is expected to have sufficient sensitivity to measure the 
X-ray spectra of sources as faint as ~ 10^^^ erg s^^ cm~^ 
in the 0.5-2 keV energy range, while the photometric limit- 
ing sensitivity is expected to be ~ 10~^® erg s~^ cm~^. In 
the hard-X band, the limiting sensitivity, both spectroscopic 

t www.rssd.esa.int/index.php?project=XEUS 
■t constellation.gsfc.nasa.gov 



and photometric, will be larger by almost an order of mag- 
nitude. At the spectroscopic flux limit of XEUS, we predict 
almost 5 x 10^ (300) AGN in the soft (hard) X-ray band, 
within a 1 deg^ f.o.v.. At such flux limits, ~ 3 x 10^ (~ 85) 
sources deg"'^ are type I objects, indicating that, because 
of obscuration, deep surveys in the soft (hard) X-ray band 
will miss nearly 90% (20%) of typo II AGN. At the photo- 
metric flux limits, we expect ~ 10® (~ 10*) sources deg~^ in 
the soft (hard) X-ray band. In this case the type II missing 
fraction is 70% (7%). 

XEUS wiW be directly probing the lower end of the mass 
function of accreting MBHs at z > 6, Mbh ~ 10®"^ Mq 
(i.e. for luminosity Lx > W^^ erg s^^ in the rest-frame 2- 
10 kcV energy band). The main contribution to the number 
counts is still expected from sources at 2 < 10, but almost 
10^ deg-^ sources, i.e., 1% of the sources, are expected to 
be detected at z > 10. We note here that, assuming an 
half energy width of the PSF of 2", XEUS wiU bo confu- 
sion limited at a sensitivity of ~ 4 x 10~^* erg s^^ cm~^, 
saturating the logN/logS at a level of ~ 2 x 10® deg"^ (Ar- 
naud et al. 2000). For bright sources, redshift determination 
should be possible from the detection of X-ray emission lines 
(e.g., 6.4 keV Ka). For fainter sources, however, redshift de- 
termination in the X-rays might bo very challenging, and 
a combination of deep X-ray observations and ultra-deep 
optical/near-IR spectroscopy is probably required (see sec- 
tion 5.3). 

In conclusion, we find that the next generation of X- 
ray missions will be able to investigate the early stages of 
MBH build-up and to provide fundamental information on 
the faint-end of the luminosity function of AGN even at very 
high redshift. 

5.2 X— ray background 

According to Moretti ct al. (2003), the intensity of the to- 
tal X-ray background (XRB) is 7.53 ± 0.35 x 10"^^ and 
2.02 ± 0.11 X 10"^^ erg s"^ cm"^ deg"^ in the 0.5-2 keV, 
and 2-10 keV energy bands, respectively. A large fraction, ~ 
94%, of the XRB in the 0.5-2 keV band has been attributed 
to sources with fluxes exceeding 2.4 x 10~^^ erg s^^ cm~'^, 
while ~ 89% of XRB in the 2-10 keV band is resolved 
into sources whose flux is > 2.1 x 10~^® erg s"'^ cm~^ 
(Moretti et al. 2003). More recently, Hickox & Marke- 
vitch (2006) estimate the unaccounted fraction of the XRB 
due to extragalactic unresolved sources as 1.77 ± 0.31 x 
10~^^ erg s~^ cm~^ deg~^ in the soft X-ray energy band 
(0.5-2 keV) and 3.4 ± 1.7 x 10"^^ erg s'^ cm'^ deg-^ in 
the hard X-ray energy band (2-8 keV). 

We compute the contribution to the residual unresolved 
XRB in the soft and hard (0.5-2 keV and 2-8 keV, re- 
spectively) energy band from the population of AGN pre- 
dicted by our model to exist at 2 > 6. Note that almost 
all these sources are below the source detection limit used 
by Moretti et al. (2003). Results are shown in Fig. 2, where 
the cumulative contribution to the XRB of sources at red- 
shift > z is considered. Different line-styles refer to differ- 
ent MBH mass ranges, as in Fig. 1. The integrate contri- 
bution to the soft XRB from AGN at z > 6 is found to be 
~ 0.4x 10~^^ erg s~^ cm~^ deg~^, corresponding to ~ 5% of 
the total, and to ~ 23% of the unaccounted fi:action. Sources 
with masses < 5 x 10® Mq give the largest contribution to 
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Figure 2. Cumulative contribution to the unresolved XRB (0.5- 
2 keV band, top panel and 2—8 keV band, bottom panel) as a 
function of redshift in erg cm~^ deg~^. Line-style as in Fig. 
1. The current observational limits (Hickox & Markevitch 2006) 
are shown as a shaded area. 



the XRB, whereas MBHs with M > 5 x 10® Mq produce a 
neghgible background. Moreover, the major contribution to 
the unresolved X-ray background is from MBHs shining at 
z > 9. In the 2-8 keV band, the XRB from unresolved, 2 > 6 
AGN is ~ 0.87 x 10~^^ erg s~^ cm~'^ deg~'^, corresponding 
to ~ 5% of the total observed hard XRB, and to ~ 25% of 
the unresolved fraction. MBHs with masses < 10® Mq give 
the main contribution also in the hard band, indicating that 
the unresolved fraction of the XRB can be used to constrain 
the faint-end of the X-ray luminosity function of AGN at 
very high redshift. 

We find that AGN at 2: > 6 contribute significantly 
to the unaccounted XRB, although their contribution is 
still well below the available constraints. Other faint un- 
resolved X-ray sources at 2 < 6 may contribute to the X-ray 
background, including galaxies, starbursts (e.g., Bauer et al. 
2004), and a population of faint AGN (e.g., Volonteri, Sal- 
vaterra & Haardt 2006; see Table 1). 

It is interesting to note that at least ~ 3 — 4% of the ob- 
served XRB (0.5-10 keV) will remain unaccounted even after 
XEUS observations will become available, as due to sources 
at 2 > 6 below the flux detection limit. Only even more sen- 
sitive X-ray observatories such as the proposed Generation- 
}3 will assess these extremely faint sources. 



5.3 JWST number counts 

Deep-field observations in the near-IR might be able, in 
principle, to detect AGN at very high redshifts. The i- 
band AGS data in the Hubble Ultra Deep Field (HUDF) 
reach sensitivity as the 28th mag (Beckwith et al. 2006). 



Figure 3. Predicted number counts of AGN at 2 > 6 in the 
1 — 3.5 fim band. Solid vertical line shows the expected JWST 
source detection limit. Lines are the same of Fig. 1. 



Our model predicts that, at this limiting magnitude, the 
sky density of AGN with 2 > 6 is ~ 0.03 arcmin"^. For 
comparison, the observed sky density of 2 > 6 galaxies in 
the HUDF is ~ 4.7 arcmin"^; Bouwens et al. 2005), corre- 
sponding to an AGN fraction of ~ 1%, so that the probabil- 
ity to detect these sources in a field as small as the HUDF 
(~ 12 arcmin^) is extremely low. Given the actual size of the 
HUDF (Bouwens et al. 2005), an area three times larger has 
to be observed at the same magnitude limit in order to detect 
one AGN at 2 > 6. In order to have a statistically significant 
number of very high redshift AGN, lower detection limits are 
needed. The James Webb Space Telescope (JWST^) IS ex- 
pected to detect sources down to a flux limit of ~ 1 njy 
(i.e., ruAB ~ 31.4), in a f.o.v. of 4 x 4 arcmin'^. We have 
computed the number counts predicted by our model for 
near-IR deep-field observations. Results are shown in Fig. 3. 
At the expected sensitivity of JWST almost 7 x 10^ deg~^ 
sources at 2 > 6 should be detected, down to a mass limit for 
the accreting MBHs of ~ 3 x 10^ Mq, almost the same mass 
scale accessed by the deepest future X-ray surveys. In spite 
of this good result, we have to note that J observations 
will be able to reveal only the un-absorbed population of 
AGN, whereas X-ray survey will probe the entire AGN pop- 
ulation. Moreover, the detection of high redshift AGN in the 
deepest JWST data could be hampered by the difficulty to 
select high-2 AGN from the more numerous population of 
high redshift galaxies. 

In conclusion. X-ray surveys are a better tool to inves- 
tigate the early stage of MBH assembly, allowing to detect 
all (type I and type II) AGN to the highest redshifts and to 
the lowest mass range. JWST deep surveys of AGN might 
be useful to investigate the optical part of the SED, and 
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Figure 4. Top panel: UV background at the Lyman Limit (912 
A) as a function of redshift in units of 10~^^ erg cm~^ Hz^^ 
sr~^. As reference value we show with the blue solid line the UV 
background as computed with the code CUBA (Haardt & Madau 
2006, in preparation) extrapolating the observed AGN luminosity 
function. Bottom panel: the corresponding HI ionization rate in 
Lines are the same of Fig. 1. 



the evolution of the type I/type II ratio at very high red- 
shifts. However, optical/near-IR follow up (with JWST and 
the forthcoming 30 meter class telescopes) of X-ray selected 
faint sources can be used to redshift determination. Un- 
obscured AGN should be recognizable both from the pres- 
ence of broad emission lines and distinctive optical-near-IR 
colours while obscured AGN should in general be identifiable 
from the detection of high- excitation optical/near-IR emis- 
sion lines. For z > 6 AGN the strongest line will be Lyman a 
at 1216 A. Further in the infrared CIV (1549 A), Hell (1640 
A), cm (1909 A) and Mgll (2800 A), will be all detectable 
by JWST, provided the gas has substantial metallicity. Of 
these lines the CIV line is probably the strongest. A combi- 
nation of JWST, 30 meter class telescopes and deep X-ray 
observations should provide the most effective identification 
strategy of AGN. 

We have to note here that MBHs shining at high red- 
shifts, similarly to miniquasars (Salvaterra, Haardt & Fer- 
rara 2005), are not expected to contribute significantly to 
the cosmic background in the near-IR, and do not represent 
a viable solution of the near-IR background excess problem 
(Salvaterra & Ferrara 2003, 2006). 

5.4 UV background and reionization 

We have also computed (see eqs. 3 and 4) the background 
intensity at the Lyman Limit (912 A) as a function of zq 
using the opacity of Fardal, Giroux & ShuU (1998; model 
Al). Results are shown in Fig. 4. As a reference, we show 
the UV background (blue solid line) as computed with an 
updated version of the code CUBA (Haardt & Madau 1996), 
extrapolating the observed AGN luminosity function and 



Figure 5. Redshift evolution of the filling factor of HH regions, 
QhiIi for different values of the clumping factor: C = 1 (dotted 
line), C = 10 (solid line), and C = 30 (dashed line). Complete 
reionization is reached when Qhii = 1- 



IGM opacity. We find a good agreement between the two 
different approaches. In the bottom panel of Fig. 4 we show 
the corresponding HI ionization rates. The UV background 
is dominated by the most massive, i.e. > 10® Mq, MBHs up 
to 2; ~ 7. At higher redshifts, these objects are rarer, and 
the background is mostly due to MBHs in the mass range 
lO'^ - lO*' Mq. At z > 10 the background intensity falls 
rapidly. 

Finally, we may ask whether high redshift AGN can 
contribute significantly to the reionization of the Universe. 
In order to answer this question, we compute the redshift 
evolution of the filling factor of HH regions as (Barkana & 
Loeb 2001), 



Qmi{z) 





dt 


/ dz' 




' z 


dz 



1 dn^, F(z',z) 

A T . ^ 1 



(5) 



where — Xhtl^ and n% are the present-day number den- 
sities of hydrogen and baryons {Xh = 0.76 is the hydrogen 
mass fraction), and dn-,/dt is the production rate of ionizing 
photons. The function F{z',z) takes into account the effect 
of recombinations. Assuming a time-independent volume- 
averaged clumping factor C, common to all HH regions, we 
can write, 



F{z',z) 



CtBriH 



C[f{z')-f{z)], 



and 



f{z) = j{i+zr + 



-13„,„3„-1 



(6) 



(7) 



where q_b ~ 2.6 x 10^ cm s~ is the hydrogen recombina- 
tion rate. 

The evolution of Qhii as function of redshift is shown 
in Fig. 5 for C = 1 (dotted line), C = 10 (solid line) and 
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Figure 6. Predicted luminosity function of quasar in the rest- 
frame liard X-ray band [2-10 keV] at z = 6. The open triangle 
shows the estimated number density of quasars in the Chandra 
Deep Field North (Barger et al. 2003b). Top panel shows the 
result for the Eddington— limited model, whereas the bottom panel 
shows the result for the rapid-growth model. 

C = 30 (dashed line). For C = 1, AGN alone are able to 
reionize completely the Universe at z ~ 7. For higher values 
of C, although the model fails to reach Qhii = 1, we find 
that the AGN contribution is not negligible and might help 
in substaining the ionization of the IGM at 2: ~ 6. 



6 COMPARISON WITH RAPID-GROWTH 
MODEL 

In this section, we discuss possible differences between the 
Eddington-limited model and a model allowing MBffs to 
accrete at super-critical rate during the early phases of 
their evolution. In Fig. 6 we plot the LF in the rest-frame 
hard X-ray band [2-10 keV] at 2 = 6 for the two mod- 
els (top panel: Eddington-limited model; bottom panel: 
rapid-growth model). Since both models reproduce the ob- 
served optical LF at 2 = 6 (Volonteri & Rees 2005, Begel- 
man et al. 2006), the X-ray LFs are also quite similar. 
For the Eddington-limited model the differential X-ray LF 
can be described by a simple power-law with power index 
j3 = —4.55, whereas for the super-critical accretion model 
/3 = —4.77. The open triangle in Fig. 6 shows The estimated 
number density of 2 ~ 5.7 AGN as obtained by Barger et 
al. (2003b) in the Chandra Deep Field North. Both models 
are in good agreement with the observed number density. 

Although the two models share similar results at 2 = 6, 
the LF at higher redshift shows significant differences. For 
example, the number density of bright quasars at 2 = 10 
(i.e. with luminosity in the rest-frame hard X-ray band in 
the range 5 x 10*^ < Lx < 3 x lO'*'^ erg s~^, correspond- 
ing to MBHs with lO'^ < Mbh < W Mq) is ~ fO"^ 
Mpc~^ for the super-critical accretion model, whereas in 



the Eddington-limited model this density decreases by al- 
most an order of magnitude. This result is confirmed by 
Fig. 7, where the solid line shows the redshift distribution of 
sources with fluxes above 10~^^ erg s~^ cm~^ (the planned 
spectroscopic flux limit of future X-ray missions) in the ob- 
served soft X-ray band for the two models. A model that 
allows an early phase of super-critical accretion will have a 
redshift distribution of sources pushed towards higher red- 
shifts with respect to the Eddington-limited model. Dotted 
lines in Fig. 7 show the redshift distribution of sources with 
X-ray fluxes above 10~^^ erg s~^ cm~^ that may also be ob- 
servable with future near-lR facilities such as JWST (i.e., at 
the source detection limit of ~ 1 njy), showing that fu- 
ture X-ray observations along with optical (near-lR) identi- 
fications by the next generation of space and 30-meter class 
telescopes might be able to discriminate the two different ac- 
cretion histories. Also shown in Fig. 7 is the redshift distribu- 
tion of sources with fluxes above 1Q~^® erg s~^ cm~^, easily 
achieved by deep Chandra observations (e.g., Alexander et 
al. 2003). Furthermore, we note here that at the planned 
soft X-ray photometric limit of XEUS, only ~ 1% of ob- 
servable sources is expected to lie at 2 > 10, if BHs growth 
is Eddington-limited, whereas this fraction increases up to 
~ 18% in case of super-critical accretion. 

ft is unlikely that sufficient 2 > 6 AGN will be identi- 
fied by the Chandra deep field observations to discriminate 
between these two models. However, we can provide indirect 
constraints by considering the contribution that each make 
to the unresolved X-ray background. We find that the con- 
tribution to the unresolved XRB increases by a factor ~ 3: 
the predicted unresolved background from sources at 2 > 6 
is ~ 1.2 X 10"^^ erg s"^ cm"^ deg"^ in the 0.5-2 keV band, 
and ~ 2.3 x 10"^^ erg s"^ cm"^ deg"^ in the 2-8 keV band. 
Results are summarized in table 1. As a net result, once that 
the contribute to the unresolved XRB of sources at 2 < 4 is 
considered, rapid-growth model saturates the observed un- 
accounted background. Note that, according to Worsley et 
al. (2005), most of the unaccounted 2-8 keV XRB actually 
lies at energies > 6 keV, strengthening our conclusions. 

We also find that the X-ray background is saturated in 
a model in which massive seeds form late, as in Koushiap- 
pas et al. (2004). In this model seed MBH form from the 
low angular momentum tail of material in halos with effi- 
cient gas cooling. In first approximation, the Koushiappas 
et al. (2004) model indicates that seed MBH form in halos 
with mass above the threshold Mh 10''Mq(1-|-2/18)"^''^, 
with a mass mscod 5 x lO*Af0(MH/lO''M0)(l + 2/18)^/^ 
The largest seed mass allows for a less efficient accre- 
tion, in particular, we have assumed that the accretion 
rate is Eddington-limited, and the maximum radiative ef- 
ficiency is 20% (Gammie et al 2004). Also in the late for- 
mation scenario, the unresolved XRB is saturated, being 
~ 1.5 X 10"^^ erg s'^ cm"^ deg"^ in the 0.5-2 keV band, 
and ~ 3.2 x 10"^^ erg s"^ cm"^ deg"^ in the 2-8 keV band. 
XRB saturation could be a problem for such models, as a 
substantial contribution to the XRB from low-redshift (i.e. 
2 < 4), unresolved faint sources (Volonteri, Salvaterra & 
Haardt 2006) and/or galaxies (Bauer et al. 2004) can not be 
excluded. 
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XRB 

0.5-2 keV 2-8 keV 


Reference 


Faint Af.;N ( 11 




\'nl.>ntrri al. 2001)" 








AGNs at z > 6 


Eddington limited 
Rapid growth 
Massive seeds 


0.4 0.9 
1.2 2.3 
1.5 3.2 


this paper 
this paper 
this paper 




Observed unresolved 


XRB 




1.77±0.31 3.4±1.7 


Hickox &6 Markevitch 2006 



Table 1. Contribution to the unresolved XRB from different sources 
in units of 10~^^ erg s~^ cm~-^ deg~^. 



« model Illb 

^ based on extrapolation down to logFi, = —18 in c.g.s units 
° based on extrapolation down to logF^, = —17 in c.g.s units 



than the HUDF area is required in order to detect just one 
of these high-z sources. 

Future X-ray missions, such as XEUS, and near-IR fa^ 
cilities such as JWST, will have the technical capabilities to 
detect accreting MBHs at z > 6 down to a mass Hmit as 
low as lO'' - 10*^ Mq. Since JWST observations will reveal 
only the un-absorbed AGN population, X-ray deep surveys 
are definitely the best tool suited for investigation of the 
early stages of MBH assembly. In particular, the next gen- 
eration of X-ray telescopes might detect as many as ~ 10^ 
sources/deg^ at their photometric sensitivity limits in the 
0.5-2 keV band. Further constraints on the population of 
high redshift AGN will be provided by future 21cm experi- 
ments like LOFAR (Rhook & Hachnolt 2006). 

We have shown that our predicted population of high 
redshift AGN would account for a significant fraction of the 
unresolved XRB (0.5-8 keV). Almost 5% of the measured 
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Figure 7. Redshift distribution of sources with fluxes in the ob- 
served soft X-ray band [0.5-2 kcV] above 10""'^'' erg s"-'^ cm~^ 
(i.e., the planned spectroscopic flux limit of future X-ray mis- 
sions, solid line), and above 10^-'"*' erg s^^ cm^^ (easily achieved 
by Chadra deep field observations, dashed line). Dotted line the 
ones with near-IR fluxes above the planned JVKSTsensitivity. Top 
panel shows the result for the Eddington-limited model, whereas 
bottom panel the one for the rapid-growth model (see Section 6) . 



7 SUMMARY AND CONCLUSIONS 

In this paper we have assessed, using Monte Carlo sim- 
ulations of DM halo merger history coupled with semi- 
analytical recipes for the assembly of MBHs within galaxy 
spheroids (Volonteri et al. 2003), the possibility of constrain- 
ing the AGN population at a > 6 with currently avail- 
able and planned space borne missions. In particular, we 
have considered ultra deep X ray and ncar-lR surveys. Wo 
claim that, among the unresolved sources in Chandra deep 
fields, few of them are AGN at « > 6, whose BH masses are 
> 10^ Mq, while in the near-IR, an are three times larger 



sources at « > 6. Moreover, in our model of 2; > 6 AGN, the 
major contribution actually comes from sources at 2; > 9, 
with fluxes < 10~^® erg s~^ cm~^, i.e. that can not be de- 
tected even by the next generation of X-ray telescopes; we 
find that at least 3-4% of the measured XRB should remain 
unresolved even after XEUS observations become available. 
These constraints become much more severe for a model in 
which super-critical accretion is allowed in the early stages 
of the MBH growth (Volonteri & Rees 2005). The model 
reproduces well the observational constraints at 2; = 6 (op- 
tical and X-ray LFs), and predicts a more extended tail 
of sources observable at 2; > 9. However, we find that in 
the rapid-growth model, the predicted XRB is three times 
higher than in the Eddington-limited accretion model, satu- 
rating the unresolved fraction of the XRB in both the 0.5 -2 
keV and 2-8 keV energy bands. Since faint sources at 2; < 4 
are expected to contribute substantially to the unaccounted 
XRB (Bauer ct al. 2004; Volonteri et al. 2006), this result 
suggests that the occurrence and effectivity of super-critical 
accretion should be investigated in much more detail. In par- 
ticular, super-Eddington accretion could be much less effi- 
cient in few-cr peaks halos, due to the gas evacuation from 
the ionizing radiation emitted by the MBH seed Pop 111 star 
progenitor (Johnson & Bromm 2006). Super-Eddington ac- 
cretion consequently is biased towards the highest density 
peaks, which experience the largest number of mergers with 
halos containing pristine gas to replenish the gas reservoir. 
Models in which seeds are much more massive than Pop 
III star remnants, as in Koushiappas et al. (2004), saturate 
the unresolved fraction of the XRB in both the 0.5-2 keV 
and 2-8 keV energy bands as well. We note here that dif- 
ferent models for MBH seed formation, although predicting 
rather large BH masses, can have a lower formation effi- 
ciency, which can ease the constraints given by the XRB 
(see Eisenstein & Loeb 1995, and Begelman, Volonteri & 
Rees 2006). Our constraints on theoretical models are con- 
servative, as we adopted the largest value of the unresolved 
XRB available in literature. 

Finally, we have computed the evolution of the UV 
background produced by the modeled population of high 
redshift AGN. Later than z ~ 7, the ionizing intensity is 
dominated by relatively massive BHs, M > 10^ Mq, while 
lighter BHs contribute mostly at earlier epochs. The UV 
background from AGN rapidly declines at z > 10. We com- 
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pute the contribution of these sources to the reionization of 
the Universe, showing that AGN alone can provide enough 
ionizing photons only in the (improbable) case of an almost 
completely homogeneous inter- galactic medium. Neverthe- 
less, for a more clumpy medium, the AGN contribution to 
the ionizing background is never negligible. 

We note that high redshift AGN can not contribute sig- 
nificantly to the near-IR background. 
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